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Abstract 

Interacting physical systems in the neighborhood of criticality (and 
massive continuum field theories) can often be characterized by just 
two physical scales: a (macroscopic) correlation length and a (micro- 
scopic) interaction range, related to the coupling and measured by the 
Ginzburg number G. A critical crossover limit can be defined when 
both scales become large while their ratio stays finite. The corre- 
sponding scaling functions are universal, and they are related to the 
standard field-theory renormalization-group functions. The critical 
crossover describes the unique flow from the Gaussian to the nonclas- 
sical fixed point. 
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Every physical situation of experimental relevance has at least two scales: 
one scale is intrinsic to the system, while the second one is related to exper- 
imental conditions. In Statistical Mechanics (SM) the correlation length £ 
is related to experimental conditions (it depends on the temperature), while 
the interaction length (Ginzburg parameter) is intrinsic. The opposite is true 
in Quantum Field Theory (QFT): here the correlation length (inverse mass 
gap) is intrinsic, while the interaction scale (inverse momentum) depends 
on the experiment. Physical predictions are functions of ratios of these two 
scales and describe the crossover from the correlation-dominated (£/G or 
p/m large) to the interaction-dominated (£/G or p/m small) regime. In a 
properly defined limit they are universal and define the unique flow between 
two different fixed points. 

In this discussion we will consider the crossover between the Gaussian 
point where mean-field predictions hold (interaction-dominated regime) to 
the standard Wilson-Fisher critical point (correlation-dominated regime). 

Massive continuum field theory is the natural setting for a description 
of this critical crossover behavior, not only in QFT, where only two scales 
characterize (super)renormalizable theories, but also in SM, where in princi- 
ple one might expect many scales (lattice spacing, geometry of interactions, 
...) to play a role, and universality may be questioned. As we will discuss, 
critical crossover scaling exists and is universal when two scales become very 
large with respect to any other (microscopic) scale. Their ratio becomes the 
(universal) control parameter of the system, whose transition from to oo 
describes the critical crossover. 

In recent years there has been extensive work H, f| H H !> H> HI 



lT| aiming at the identification of the correct, theoretical and experimental, 
definition of the critical crossover limit. We give here a sketch of the argument 
for (i-dimensional iV-component vector spin models, but the notion may easily 
be extended to many other physical systems. 

Let us start with the standard Landau- Ginzburg Hamiltonian on a d- 
dimensional lattice, 

Tt = ^ \J{xj - xj) [4>{xj) - (pjxj)} 2 



where </>(xi) are iV-dimensional vectors. We will first consider the short- 
range case in which J(x) is the standard nearest-neighbour coupling. For 
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this model the interaction scale is controlled by the coupling u and the rel- 
evant parameters are the (thermal) Ginzburg number G and its magnetic 
counterpart Gh defined by: 

G = U 2 '^- d \ G h = ^+2)/[2(4-d)]_ (2) 

Under a renormalizat ion-group (RG) transformation G scales like the (re- 
duced) temperature, while Gh scales as the magnetic field. For t = r—r c <C G 
and h ^ Gh one observes the standard critical behaviour, while in the oppo- 
site case the behaviour is classical. The critical crossover limit corresponds 
to considering t, h,u — > keeping t = t/G and h = h/Gh fixed. This limit 
is universal, i.e. independent of the detailed structure of the model: for 
Hamiltonians ([!]) the same behaviour is obtained as long as the interaction is 
short-ranged, i.e. for any J(x) such that J2 x x 2 J(x) < +00. The crossover 
functions can be computed in the standard continuum 4 theory 0, |], ||. 
A dimensional analysis shows that (using the subtracted bare mass and re- 
moving the cutoff) finite results can be obtained directly in terms of the 
dimensionless variable u/t 2 ~ d ^ 2 = t d ^ 2 ~ 2 , and no further limiting procedure 
is required. It is important to observe that the critical crossover functions 
are related to the standard continuum RG functions if one expresses them 
in terms of the zero-momentum four-point renormalized coupling g. The 
crossover functions are well studied H |], || in the fixed- dimension expansion 
when d = 3. 

Let us now consider the long-range case. We assume that J(x) has the 
following form 

where D is a lattice domain characterized by some scale R. Explicitly we 
define R and the corresponding domain volume V# by 

vk-EL * = m£*- (4) 

x<=D n x£D 

The shape of D is irrelevant for our purposes as long as Vr ~ R d for R — > 00. 
The constant J defines the normalization of the fields. In our discussion it 
is useful to assume a long-range normalization, i.e. J = 1/Vr, since with 
this choice the limit R — > 00 is well-defined. Notice the this is not the 
normalization that is commonly used discussing short-range models. Indeed, 
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in the latter case, one defines J = R~ 2 /Vr, so that the propagator behaves 
as k 2 for k — > 0. To understand the connection between the theory with 
long-range interactions and the short-range model let us perform an RG 
transformation ||. Define new ("blocked") coordinates yi = Xi/R and rescale 
the fields according to 



ft Vi ) = R d/2 (t>{Ryi), h{ Vi ) = R d/2 h(R yi ). 
The rescaled Hamiltonian becomes 



(5) 



2 



'■J 



(6) 



where now the coupling J(x) is of short-range type, i.e. independent of 
R. Being short-ranged, we can apply the previous arguments and define 
Ginzburg parameters: 

2/(d-A) 



G = (uR~ d ) 
G h = R- d ' 2 (uR- d ) 

= u (d+2)/[2(d-4)] R -3d/(4-d) 



u 2/(d-4) R -2d/(4- 
d , (d+2)/[2(d-4)] 



d) 



(7) 



(8) 



Therefore, in the long-range model, the critical crossover limit can be defined 
as R — > oo, t, h — > 0, with t = t/G, h = t/G^ fixed. The variables that 
are kept fixed are the same, but a different mechanism is responsible for the 
change of the Ginzburg parameters: in short-range models we vary u keeping 
the range R fixed and finite, while here we keep the interaction strength u 
fixed and vary the range R. The important consequence of the argument 
presented above is that the critical crossover functions defined using the 
long-range Hamiltonian and the previous limiting procedure agree with those 
computed in the short-range model, apart from trivial rescalings. 

Let us give a few examples. Let us introduce magnetic susceptibility, 
correlation length and magnetization in the usual way: 



X 

v X 

M = (0 O >. 



(9) 

(10) 
(11) 
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Then one can show that in the limit t — > 0, G, Gh — > with t and h fixed 
the following rescaled quantities have a finite limit: 

X = X G^F x (t,h), ^ (12) 
e = R- 2 eG^F e (t,h), (13) 
M = MG/G h — > F M (t,Ji), (14) 

where F x (t,h), F^(t,h), and F M (t,h) are universal apart from a rescaling 
of t and h and of the functions themselves. Comparison with experimental 
data is usually performed introducing effective exponents. For instance, for 
h = 0, we can define 

7e ff (t) = -filogi^ftO), (15) 
at 

u eS (t) = ---^logi^O). (16) 
* at 

These functions can be related to the standard RG functions if one expresses 
them in terms of the zero- momentum four-point coupling g. In the high- 
temperature phase one finds for instance (cf. Ref. H) 

lesjg) = 7(5) ( 17 ) 

^eff(^) v(gY 

v{g)m d -^- = i{g)-iM, (18) 

where 7(^7), ^(g), and (3(g) are the standard RG functions. These effective 
exponents interpolate between the classical and the non-classical value. As 
an example, in Fig. [I], we report the graph of 7 e ff(£) in the high- and low- 
temperature phase for the Ising model (the computation has been done using 
the results of Refs. 0, f|, |fl). These curves are in good agreement with 
numerical results for the long-range Ising model [|12||, even for very small 



values of R, i.e. for interactions extending over a few lattice spacings. 

The ideas we have presented here can be explicitly checked in the large- 
iV limit. All the crossover functions can be computed in the whole (t, h)- 
plane for 2 < d < 4 and universality (model independence) can be explicitly 
checked. For instance for the effective exponents defined in Eqs. (|Uj) and 
(P~6| ) one finds in three dimensions 

7eff (t) = 2v cS (t) = 1 + (1 + aft- 1 ' 2 , (19) 
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where c 7 is a non- universal constant. 

One can also study the corrections to the leading universal behaviour. If 
R is chosen as in Eq. (|J), one verifies that the corrections to the univer- 
sal crossover functions scale as 1/R d (for generic choices of scale one would 
observe instead l/i?-corrections). 

The discussion of these non-universal effects can be extended to all values 
of iV considering a perturbative expansion around the mean-field solution. 
Consider for instance the long-range Ising model 

H= -^-^J{xi-Xj)s{xi)-s{xj), (20) 

where J(x) is defined in Eq. (^). For N = 1 and for a particular choice of 
D this is the model that has been studied in Refs. || ||, pL2|| . Computing the 
corrections to mean field allows us to determine the corrections to f3 c (R) for 
R — > oo for the Hamiltonian (|20|). One finds for d > 2 

/r\ d h 1 
W?W) (21) 

with corrections of order R~ 2d ( multiplicative logarithms appear for some 
special values of d, for instance for d = 3). Here 

n ^ = 7T E i 1 - ^ > ( 22 ) 

and the integral is extended over the first Brillouin zone. Expanding the 
integral in powers of R one finds 

M R ) = 1 + ^ + • • • (23) 

where a depends on the precise definition of the domain D. In two dimensions 
there are logarithmic corrections and one finds 

Pe{R) = 1 + log^ 2 + 0(R- 2 ). (24) 

By considering the mean-field limit one can also relate the non-universal 
constants that appear in the definition of the crossover scaling functions. To 
give an example, consider x for h = 0. The function F x (t, 0) can be computed 
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in perturbation theory, cf. Ref. [|J, obtaining a function F^(t). On the other 
hand, if one considers the model ( |20| ) one obtains a different F^(t). As we 
explained above we should have 



F x(t)= "xW- ( 25 ) 



The analysis of the mean-field limit provides exact expressions for a x , which 
depends on the observable, and b, which depends only on the model. 

Finally we want to discuss the crossover scaling limit for models that 
have (3 C = +00. This is the case of the two-dimensional N- vector model with 
N > 3. For these theories define 



and consider the limit R — > 00 with t fixed. One finds that the limits 
defined in Eqs. (|12"D and (|T3"|) still exist and define crossover functions of 
t. For t — ► +00 these functions show mean-field behaviour, while standard 
asymptotic scaling is observed for t — > — 00. Notice that one can use t as 
a variable instead of t also when /3 C is finite. In this case, however, nothing 
new is obtained, since t — t is simply a constant for R — > 00. 

The model with Hamiltonian (|20|) can be studied in the limit iV — > 0. In 
this case it can be rewritten in terms of self- avoiding walks (SAWs) [13j with 
long-range jumps. To be explicit, we define an n-step SAW with range R as a 
sequence of lattice points {u , • • • , u n } with uj = (0, 0, 0) and ujj + \ e Dn(uj), 
such that uJi 7^ cjj for all i 7^ j. Then, if c n ^(x) is the number of n-step SAWs 
with range R going from to x, we indicate with c nj # the total number of 
n-step walks and with E\ R the mean square end-to-end distance. They are 
defined as: 




(26) 




(27) 



X 



El 




(28) 



One can then prove that 



00 




(29) 



n=0 




n=0 



(30) 
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where (3 = (3/Vr and x an d £ 2 are defined in the model (PH|). 

The crossover limit is trivially defined remembering that n is the dual 
variable (in the sense of Laplace transforms) of t. Therefore we should study 
the limit n — > oo, R — > oo with n = nR~ 2d ^ 4 ~ d ^ fixed. From Eqs. (|T2"D and 
(|j~3| ) we obtain that the following limits exist: 

c„,fi= Cn >R (3 c (R) n -> <? c (n), (31) 
^ = ^, R i?- 8/{4 - d) - ^(n), (32) 

where the functions <7 c (ft<) and g E {n) are related by a Laplace transform to 
F x (t,0) and Fp(b,0). Explicitly 

/"OO 

F x (t,0) = / dug c {u)e- ut , (33) 
Jo 

1 f°° 

F e {t,0)F x {t,0) = Ydj Q du 9c(u)g E (u)e- ut . 



(34) 

Using perturbation theory it is possible to derive predictions for E\ and c n . 
For E\ we can write 

gE,pr(n) = a E nh E (z), (35) 

where z = (n/l) 1 ^ 2 . The function h E (z) has been computed in perturba- 
tion theory to six- loop order []14| . Resumming the series with a Borel-Leroy 
transform one finds that a very good approximation is provided by [JT5J 

^(2) = (1 + 7.6118* + 12.05135^ 2 ) - 175166 . (36) 

Comparison with a detailed Monte Carlo simulation for the Domb- Joyce 



model indicates |1| that this simple expression differs from the exact result 
by less than 0.02% for z < 2. 

The constants a E and I appearing in Eq. ( |35D are non universal. For our 
specific model they are given by 

a E = 6, I = (4tt) 3 . (37) 

We have performed flnj an extensive simulation of this model of long-range 
SAWs generating walks of length up to iV pa 7 ■ 10 4 . The domain D was 
chosen as follows: 

D= Ix-.V] \ Xi \ <p \ . (38) 
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In the simulation we varied p between 2 and 12. Let us describe the results 
for the end-to-end distance. Analogous results can be obtained for c n ^. 

In Fig. H we report our results for E% R together with the perturbative 
prediction gE,pr{p) defined in Eqs. (|35],[3B]). The agreement is very good al- 
though one can see clearly the presence of corrections to scaling. In order to 
see better the discrepancies between the numerical data and the theoretical 
prediction we report in Fig. |3| the ratio E% R /gE,pr(n)- In this plot the cor- 
rections to scaling are clearly visible: points with different R fall on different 
curves that converge to 1 as expected. For p = 12, the deviations are less 
than 0.2%. It is interesting to observe that the corrections change sign with 
n. For n < 2 ■ 10 2 , the corrections are negative, while in the opposite case 
they are positive. 

The corrections to scaling are expected^ to scale as R~ d . To check this 
behaviour let us consider 



A 



E;n,R 



El R /9E,PT{n)-l R 3 - (39) 



The plot of A.E;n,R is reported in Fig. f|. A good scaling behaviour is observed 
confirming the theoretical prediction for the corrections. This nice scaling 
indicates also that the approximation (^) can be considered practically exact 
at our level of precision. 

We have also defined an effective exponent z/ eff 

— (40) 

It is reported in Fig. [|. It shows the expected crossover behaviour between 
the mean-field value v = 1/2 and the self-avoiding walk value v = 0.58758(7) 

Pi- 
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Figure 1: Effective susceptibility exponent as a function of t for the high- 
(7eff) and low- (7^) temperature phase of the three-dimensional Ising model. 
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Figure 3: Results for E 2 R j g E ^ PT {n). We use the same symbols as in Fig. 
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Figure 4: Results for lS. E . n)R . We use the same symbols as in Fig. ^|. 
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Figure 5: Results for z/ eff . We use the same symbols as in Fig. §. The dashed 
line is the self-avoiding walk value v = 0.58758(7). The solid line is the 
theoretical prediction. 



15 



